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Outline 


• Background and Motivation: 

• Cryogenic fluid behavior in flight conditions 

• Long-term, in-space storage of cryogenic propellants for future exploration missions 

• Problem Setup: 

• Geometry & grid 

• Fluent settings 

• Fluid properties (N 2 ): temperature, pressure dependent? 

• UDF for condensation / evaporation 

• UDF for non-inertial reference frame 

• Time-dependent acceleration and Bond number 

• Boundary conditions 

• Procedure for initial conditions 

• Comparison with experimental data: 

• Initial thermal profile 

• Visual comparison with high-speed movie 

• Pressure data: balance of evaporation and condensation 

• Net heat transfer/boiling heat transfer rate 

• Comparison with temperature sensor data 
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Background & Motivation 




Background: 

• LN2 tank in 2010 low-g parabolic aircraft campaign 

• Significant condensation, evaporation, & boiling 

• Simulation compared to one low-g parabola 

Motivation: 

• Cryogenic fluid behavior in flight conditions 

• Long-term, in-space storage of cryogenic propellants for 
future exploration missions 
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Geometry & Grid 




3-Dimensional grid, 360 degree sector 
Fluid Grid: 569,110 Cells 

• In interior, uniform, structured grid 

• 1 mm resolution 
Solid Grid: 685,858 Cells 

• Unstructured grid 

• Variable resolution 

Thermal isolation at joint, sealing gasket 
Post mounted temperature sensors, not 
simulated 

Currently, no refinement 
Partitioned for 16 or 32 processors 

Fluid tank dimensions: 

• Radius: ~3 cm 

• Height : ~10 cm 
Slosh frequency: 

• Observed ~4 Hz 

• Calculated 5.0 Hz 
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Lid Surface 


Geometry & Grid 
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Lid Surface 


Geometry & Grid 
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Fluent Setup for Simulations 

• Simulations performed using ANSYS Fluent version 13; 

• 3D grid of fluid and solid regions 

• Mass, momentum, energy, turbulence PDEs 

• Compressible, ideal gas; Boussinesq liquid 

• Fluid properties of nitrogen for fluid viscosity, thermal conductivity, specific heat, latent heat 
of vaporization, surface tension at 77.244K, 1 bar from webbook.nist.gov/chemistry/fluid 

• Solid: temperature dependent density, specific heat, and thermal conductivity from CNES for 
inox (stainless steel), aluminum, and sapphire 

• Volume of Fluid (VOF) for 2-phase flow 

• k-co SST turbulence model of Menter et. al. (turbulent damping = 10) 

• UDF for mass transfer at liquid/ullage interface, boiling, and gas phase condensation 

• UDF for non-inertial acceleration, 

• Boundary conditions on later slide, 

• Second order upwind scheme was used for discretization of the mass, momentum, energy, 
and turbulence, (cell values) 

• PISO scheme was used for the pressure-velocity coupling (cell values) 

• Least Squares Cell Based scheme was used for the gradient calculations (face values) 

• PRESTO! scheme was used for the pressure interpolation (face values) 

• First order implicit time discretization was used, also bounded second order implicit 

• Time step = le-04 seconds 

• Extensive instrumentation of the simulation 
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Temperature & Pressure Dependence of Fluid Properties 


Initially used constant properties for p, k, C d , at 77.244K, 1 bar, but ~10% error 


p' 

Range of temperature and pressure: 70-110 K; 

Well represented by simple polynomial in temperature: 

• Viscosity p, (vapor & liquid) 

• Thermal conductivity k (vapor & liquid) 

• Surface tension, y 

• C p liquid 

Varies with both temperature and pressure: 

• C p gas 

• Heat of vaporization / condensation 

• Z varies ~8% (Ideal gas assumption) 
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Evaporation/Condensation UDF 


Mass transfer and heat of vaporization/condensation 


• Based on Hertz-Knudsen-Schrage equation: ^ 
• Derived from Maxwell dist'n 


2 - rw I 2nR u 


2 \MW va p / Psat(Tliq ) 


kg/s-m 2 Evaporation is +ve 


• Assume: ff cond — &evap > 



• Constant(o) * (P sat (T) - P)/sqrt(T)/ length_scale for local P, T 

• Enforces saturation conditions on interface 

• UDF Define_Adjust() calculates mass transfer; UDF Define_Mass_Transfer() applies 

• Requires kg/s-m 3 , hence length_scale = sqrt(l/ 1 grad cl 2 ), c is VOF fraction 

• Different situations, different accommodation coefficients: 

• Interface condensation, a = l.OxlO' 4 is 'best' fit 

• Interface evaporation, o = l.OxlO' 4 is 'best' fit 

• Boiling (liquid phase evap), a = 5.0xl0' 3 is 'plausible' fit 

• Gas phase condensation, a = l.OxlO' 4 is used 

• Boiling— vaporization away from a liquid/vapor interface: 

• Superheat criteria in each cell: T max -T sat (P) > 5 K T max is max in cell (walls too) 

• Dry boiling cut-off 

• P sat (T) by curve fit from Reynolds, Thermodynamic Properties in SI 

• T sat (P) curve fit to NIST data 
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Non-lnertial Reference Frame UDF 



• Non-inertial reference frame accounts for: 

• Linear acceleration of aircraft, a 

• Angular rotation, w, and angular acceleration, a (not present) 
a cg + 2£0Xt^+6JXf + 6JX(6JXf) 

• In general, piece-wise linear fit to a, w, v r samples 






Here, acceleration, a, sampled at 2 Hz, 10 Hz 

• Two components: a x , a z ; assume a y = zero 

• Piece-wise linear fit to supplied a x , a z 

• Initial conditions: steady acceleration 

a x =-16.5 m/s 2 , a y =0, a z =-1.93 m/s 2 

UDF Define_Source adds terms to RHS of 
momentum equations as p a, (kg m 2 /s 2 ), and 
RHS of energy eqn. as p a*v, (kg m 3 /s 3 ) 

Fluent has trouble with microgravity 

• Issues at 0( a/g ) ~ 10' 6 

• No issues at 0( a/g ) ~ 10" 5 
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Simulation Conditions: Bond Number 




Bond number: 

• Range [ 0.3, 6. ] 

• Mean 2. 

• 0(1) for 10 s 

• Surface tension forces/ 
body forces 

• Eotvos number 


• Acceleration due to gravity at 10 Hz supplied by CNES 

• 2 components: g x , g y 
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Boundary Conditions 



Inox 

(Stainless 

Steel) 


Sapphire 


Aluminum 



Heat fluxes total ~4W 

Assumed constant on surface of each part 

Due to radiation & conduction 

Liquid-to-vapor contact angle: 5 degrees 

Initial constant acceleration: 

• a x = -16.4993 m/s 2 ; a x /g = -1.682 

• a y = 0. 

• a z = -1.9325 m/s 2 ; a z /g = -0.197 
Initial interface 

• Position: ~60 mm from bottom 

• Angle: from initial acceleration 
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Initial Thermal Conditions 





Initial conditions by transient fluid- 
thermal simulation 

• Constant gravity 

• 90 s with time step of 4.0xl0' 4 s 
Thermal isolation at joints 
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Temperature & Heat Flux: What & Why? 



Relative Thermal Conductivity Initial Heat Flow Heat Flow in Re-Orientation 

Heat Build-Up 




Initial 

Temperature 

Distribution 



Vapor and sealing gasket create insulation (high temperature gradient) 
Inox lid heats up in high-g interval with vapor at top 
With low-g, re-orientation, liquid impinges on hot lid, and boils 
Heat Transfer: surface boiling, departure of bubble, condensation 
Heat transferred into well-mixed liquid with high heat storage capacity 
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Section 2: Comparison with Experimental Data 



• Initial temperature profile 

• Visual comparison with high-speed movie 

• Pressure data: balance of evaporation and condensation 

• Net heat transfer/boiling heat transfer rate 

• Comparison with temperature sensor data 
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Initial Thermal Conditions 





Simulation 

Temperature, 

T-T 0 

(K) 
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Difference 
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Discrepancy near tl2b in high temperature gradient: gasket modeling? 
Discrepancy near tl2l and lower lid: specified heat fluxes? 
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Initial Re-orientation of Surface 
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Heavy Boiling Phase with 
Condensation and Transit 
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Pressure Evolution 
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Internal Energy: A Measure of Heat Transfer 



Is Boiling Heat Transfer Rate Correct? 
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Results: Wall Temperatures 





o 

90 92 94 96 98 100 102 104 106 108 110 112 114 116 118 120 

Time (s) 



www.nasa.gov 22 


National Aeronautics and Space Administration 


Results: Temperature Sensor T12G 
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Fluid Temperature on Midplane 
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Results: Temperature Sensor T12A 
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Temperature Sensors: Dry-to-Wet, Wet-to-Dry 




~0.1 s ~0.1 s 0.3-0. 5 s 
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Diode sensors time 
constant: t = 0.1 s 

95% in 3 time constants, t 

Hot gas exposure duration is 
0.3 -0.5 s (one-way) 

Wet-to-dry transition 
includes a liquid film that 
must vaporize, before gas 

Wet-to-dry time delays 
observed experimentally 
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Sensors: 
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Wet-to-Dry With Drainage 



Temperature discrepancy between T12a sensor (top) and simulation 

After low gravity phase (final re-orientation), as lid should be heating 
vapor to create a stable thermal stratification 

Experimental geometry is different: fill line and valve 
Wet-to-dry transition complicated by drainage? 

Drainage of liquid visible, in experiment, 2 s after final re-orientation 
Some simulations shows waves in thermal stratification, others do not 
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Summary 

• Presentation: setup of Fluent and comparison with experimental results 

• Initial temperature profile 

• Visual comparison with high-speed movie 

• Pressure data 

• Net heat transfer/boiling heat transfer rate 

• Comparison with temperature sensor data 

• Generally, good agreement with experimental data 

• Evidence for low sensitivity of wet-to-dry temperature sensors 

• Limitations of boiling model 

• Limitation in prediction of condensation / evaporation 

• Future work: 

• Further analysis of thermal layers near fluid/vapor interface 

• Grid resolution studies 
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